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ABSTRACT 

We show how the recent discovery of a likely close white dwarf companion to the well known star 
Rcgulus, one of the brightest stars in the sky, leads to considerable insight into the prior evolutionary 
history of this star, including the cause of its current rapid rotation. We infer a relatively narrow range 
for the initial masses of the progenitor system: M\q — 2.3 ± 0.2 M & and M20 = 1.7 ± 0.2 M 0I where 
M 10 and M 2 o are the initial masses of the progenitors of the white dwarf and Regulus, respectively. 
In this scenario, the age of the Regulus system would exceed 1 Gyr. We also show that Regulus, with 
a current orbital period of 40 days, has an interesting future ahead of it. This includes (i) a common 
envelope phase, and, quite possibly, (ii) an sdB phase, followed by (iii) an AM CVn phase with orbital 
periods < 1 hr. Binary evolution calculations are presented in support of this scenario. We also discuss 
alternative possibilities, emphasizing the present uncertainties in binary evolution theory. Thus, this 
one particular star system illustrates many different aspects of binary stellar evolution. 
Subject headings: binaries: spectroscopic - stars: early-type - stars: variables - stars: dwarf novae - 
stars: individual (Regulus, a Leo) 



1. INTRODUCTION 



1.1. 



Properties of the Regulus System 

Regulus (a Leo; HD 87901) was recently discovered to 
be a spectroscopic binary which has a likely white dwarf 
companion of mass M w <j — 0.3 Mp, and an orbital period 
of 

-Porb — 40.11 d (|Gies et al.|[2008[ hereafter "the inner 
binary system"). Regulus itself has an inferred mass of 
~3.4±0.2 Mq. This very bright star has been known for 
many years to have at least two other c ompanions (BC) 
which together form a binary system (McAlister et al. 
2005 1. This BC subsystem is located 177" from Regulus, 



too great a distance to have ever directly interacted with 
the star. The B component (a Leo B; HD 87884) is a 
~0.8M© star of spectral type K2V; the C component 
is a very faint M4V star with a mass of ~0.2Af Q . The 
Washington Double Star Catalog (Mason et al. 2001) 
lists a D component of the system, also having a common 
proper motion with the system and a separation of 217" 
from component A. 

Located at distance of 24.3±0.2 pc, Regulus is a star of 
spectral type B7V with an apparent magnitude of 1.36, 
making it one of the brightest stars in the sky. The star 
rotates very rapidly, with a rotational period of 15.9 h, 
const ituting ~86% of its break-up speed ( |McAlister et al. 
2005 1. Its rapid rotation causes it to be highly flattened, 
with an equatorial diameter 32% greater tha n its polar 
diameter (4 .16 Rq and 3.14 R@, respectively; McAlister 
eFaIp005| |. 



There is som e deba te regarding the age of Regulus. 
Gerbaldi et al.| ( |2001| ) estimate it to be ~150Myr old; 
however, a comparison of its age with that of compo- 
nent B, assumed to be coeval, results in a 100 Myr dis- 
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crepancy. It is worth noting that the age estimate for 
Rcgulus is based largely on the effective temperature of 
the star. However, the rapid rotation causes a d ifference 
of ^5000 K bet ween the poles and the equator (McAlis- 
ter et al. 2005), making this assessment somewhat sus- 
Moreover, this estimate assumes that Regulus, 
3.4 M©, has evolved as a single star. As we 



pect 
with M 

will show in this paper, the original mass of Regulus was 
almost certainly much lower, and Rcgulus has attained 
its present mass by mass transfer from its close compan- 
ion. If this scenario is correct, it is difficult to escape the 
conclusion that the entire Regulus system has a probable 
age of > 1 Gyr. 

In this paper we use the fact that Regulus has a white 
dwarf companion of close to 0.3 Mq in a 40-day orbit to 
reconstruct an approximate past history for the system 
(§2). Based on the current state of the Regulus inner 
binary, we discuss various possible future channels for the 
system (§3) with one of the most interesting possibilities 
that it will become an AM CVn binary with a period as 
short as a few minutes. 

2. PRIOR EVOLUTION OF THE INNER REGULUS BINARY 



2.1. R 



orb 



M w( j Relation 



It has been known for a long time that accretion from 
a companion star can result in very high rotation rates 
for the accretor. However, in the case of Regulus, no 
close com panion star was k nown to exist until the dis- 
covery by |Gies et al. (2008 1 that Regulus has a co mpan- 
ion in a 40-d ay orbit that is likely a white dwarf. Gies 
et al. (20081 suggested that matter accreted by Regu- 
lus from the envelope of the progenitor of the putative 
white dwarf was responsible for the high rotation rate. 
Since only a mass function is measured for the Regulus 
inner binary, the mass of the u nseen companion star is 



formally limited to M > 0.3 Mq. |Gies et al.| ( |2008[ ) argue 
against a neutron star companion on the grounds that a 
~1.4M companion would require an orbital inclination 
of i < 15° which has only a small a priori probability of 
~3%. No formal limit on the orbital eccentricity is given 



2 



by Gies et al. (20081, but an inspection of their Fig. 1 
suggests that e % 0.05. This makes it even more unlikely 
that the unseen companion is a neutron star considering 
that the loss of even a small amount of matter during 
the supernova explosion that gave birth to the neutron 
star would undoubtedly have led to a sizable eccentricity. 
We present another argument below that the mass of the 
unseen companion is almost certainly near 0.3 Mq. 

Assuming that the unseen companion to Regulus is a 
white dwarf, we can show why the 40-day orbit is just 
what is expected on theoretical grounds. For stars with 
an initial mass < 2.5 M there is a nearly unique rela- 
tion between its properties on the giant branch and its 
He core mass. Since, for Roche-lobe filling donor stars 
that are lower in mass than the accretor (or even slightly 
higher), the orbital period depends only on the mass and 
radius of the donor star, the orbital period at the end 
of mass transfer from a giant star depends only on the 
core mass of the donor star. This period— white dwarf 
mass relation has been studied in detail by Rappaport 
et al. (1995). They found the following approximate an- 
alytic expression: 



P orb ~ 1.3 x 10 5 M£ d 25 (l + 4M^ d ) 



3/2 



(1) 



where the orbital period is expressed in days and the 
white dwarf mass is in units of Mq. Rappaport et 
al. (1995) estimated a theoretical uncertainty in the 
white dwarf mass of ~ ±18% that takes into account un- 
certainties in the chemical composition, the initial mass 
of the parent star, and the mixing length parameter. If 
we solve this expression for M wc j with P or b = 40.1 d, we 
find 

M wd ~ 0.28 ± 0.05 M G (2) 

(full uncertainty). This is entirely consistent with the 
measured value of M w d > 0.3 Mq. 

If we assume that the rotation axis of Regulus, which 
is measured to be > 75° with respect to the line of sight 
(McAlister et al. 2005), coincides with the normal to the 
orbital plane, then the orbital inclination angle of the 
Regulus inner binary should be i > 75°. That this is the 
case follows from the assumption that the matter trans- 
ferred from the progenitor of the white dwarf compan- 
ion was responsible for spinning up Regulus (see [2.4 1. 
Therefore, the angular momentum vector of the orbit and 
that of the rapidly rotating Regulus should coincide. For 
i > 75°, the mass function for Regulus (Gies et al. 2008) 
yields amass for the white dwarf of M w d = 0.302±0.017. 
Again, this is highly consistent with the theoretical value 
expected for a white dwarf remnant in a 40-day circular 
orbit (see eq. 2; Rappaport et al. 1995). Finally, we note 
that, according to the theoretical scenario for forming 
the Regulus inner binary, the orbital eccentricity is ex- 
pected to be e < 10~ 4 (see Rappaport et al. 1995, and 
references therein) - which can possibly be falsified in 
the future. 

2.2. Constraints on the Primordial Binary 

In the evolutionary scenario for the current Regulus in- 
ner binary, the progenitor of the white dwarf is the more 
massive star, while the secondary - the progenitor of the 
current Regulus - is somewhat less massive. The initial 
orbital period must be less than 40 days in order for the 
He core of the primary not to exceed ~0.3 Mq (see eq. 1). 
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Fig. 1.— Contours of (3, P orb , t K h,i/t K h,2, and r nuCi i/T nuc , 2 in 
the Mio — M20 plane of the progenitor of the Regulus inner binary. 
All physically realizable points in the Mio — M20 plane must lie 
above the heavy solid line with negative slope —1 (common to 
all panels) in order to yield a current binary mass of 3.7 ± 0.2 Mq . 
Similarly, all physically meaningful points must lie below the heavy 
line of equal masses with slope +1. 



When the primary overfills its Roche lobe, its envelope is 
transferred stably, though not necessarily conservatively, 
to the secondary. At first, the orbit shrinks due to the 
transfer of mass from the more, to the less, massive star. 
Once the mass ratio reaches unity, the orbit will start to 
expand. The ratio of the final-to-initial orbital period is 
given by 
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where C x = 3a(l - /3) - 3 and C 2 = -3a(l - /3)//3 - 3 
(Podsiadlowski et al. 1992). The parameter j3 is the frac- 
tion of the transferred matter retained by the secondary, 
and a is the specific angular momentum in units of the 
specific orbital angular momentum of the system carried 
away by any matter lost from the system. The notation is 
that "1", "2", and "b" stand for the primary, secondary, 
and the binary, respectively, while "f" and "i" indicate 
final and initial states, respectively. The derivation of 
eq. (3) involves the assumption that the parameters a 
and (3 are constant throughout the mass-transfer phase; 
this approximation seems justified in light of the other 
theoretical uncertainties for this phase of the evolution. 

In order to study the range of possible progenitor stars 
for the Regulus inner binary system, we considered pri- 
mary stars of mass M 10 < 3M Q , and secondaries of any 
mass less than the primary mass. For primaries with 
Mio It 2.5 — 3 Mq the core mass that develops would 
be M core > 0.5 Mq and would not leave a remnant con- 
sistent with the 0.3 Mq white dwarf that is inferred to 
be orbiting Regulus. For each point in the Mio — Mao 
plane we estimated the ratio of thermal timescales, tkh, 
and nuclear timescales, r nuc , for the primary and sec- 
ondary. Contours of constant timescale ratios are shown 
in the bottom two panels (c and d) of Fig. 1. We take 
tkh ^ M 2 /(RL) and r nuc cx M/L, with stellar luminos- 
ity and radius scaling like M 3 7 and M 8 , respectively. 
We then find tkh oc M~ 



-2.5 



-2.7 



Evolution of Rcgulus 



3 




2.4 
M7M, 



Fig. 2. — Allowed regions in the Afio — M20 plane for the pro- 
genitor of the Regulus inner binary (see Fig. 1). The cutoff for 
Afio > 2.5 Mq indicates the highest mass primary that could yield 
a white dwarf mass as low as 0.3 Mq. The black square indicates 
the initial masses used in the numerical evolution model that pro- 
duced the best match to the current Regulus (inner binary) system. 

In the top left panel of Fig. 1 (panel a), we show con- 
tours of constant /3, the fraction of transferred mass re- 
tained by the secondary, in the M\q — M20 plane. These 
values are inferred by comparing the initial mass of the 
primordial binary with the current observed mass of 
M b ,f ~ 3.4 + 0.3 Mq = 3.7 M Q . Finally the top right 
panel (b) in Fig. 1 shows contours of constant initial or- 
bital period that would lead to the observed current value 
of forb = 40 days. To compute P or b we utilized eq. (3) 
with P taken from Fig. la, and a, the specific angular 
momentum associated with mass loss, taken to be a con- 
stant throughout the mass loss proces s an d to have a 
typical fiducial value of unity (but see [2.5 1. 

If we require that the ratio of thermal timescales for 
the progenitor stars, Tkh,i/tkh,2 be not too different 
from unity, so that the secondary can retain a sizable 
fraction of the transferred matter, we can restrict the 
allowed range of progenitor masses to lie above a par- 
ticular contour in Fig. lc. We somewhat arbitrarily 
choose tkh.i/tkh,2 ^ 0.4 in order that the thermal 
timescales not be too disparate. Similarly, we require 
that the ratio of nuclear timescales for the progenitor 
stars, T nuCj i/r nuc 2 be sufficiently different from unity, so 
that the current Regulus has not already evolved up the 
giant branch. Since, as we show in §2.5, the initial sec- 
ondary star will spend ~200 Myr with sufficient accreted 
mass to exceed the mass of the original primary, and 
subsequently, another ~50 Myr with much more mass 
than the primary, as it approaches its current mass of 
3.4 Mq, the original mass of the secondary cannot be 
too close to that of the primary. In order for Regulus 
not to be more evolved than its current state, we require 
that T nuCj i/r nuc ,2 0.6, and we adopt this as a rough 
upper boundary in the bottom right panel of Fig. 1. 

2.3. Parameters of the Primordial Binary 

Combining these constraints yields the allowed range 
of M 10 and M 2 q (see Fig. 2) . The most likely primordial 
masses are M w ^ 2.3±O.2M and M 20 ^ 1.7±O.2M . 
The initial orbital period could have been anywhere in 



the range of P OI b ~ 1 — 15 days. Calculations of grids, 
and even entire populations, of such mass transfer bi- 
naries are given, e.g., in Nelson & Eggleton (2001) and 
Willems & Kolb (2004). These authors considered only 
conservative mass transfer, but their detailed evolution- 
ary calculations arc nonetheless quite instructive. 



2.4. The Origin of the Rapid Rotation 

It has long been argued that one way to produce a very 
rapidly rotating star (e.g., a Be star) is by the accretion 
of mass and angular momentum from a companion star 
(e.g., Pols et al. 1991), similar to the origin of the Be phe- 
nomenon in Be-/X-ray binaries (Rappaport & van den 
Heuvel 1982). Generally, a star has to accrete < 10% of 
its initial mass to be spun up to near critical rotation 4 
(Packet 1981). But, if the companion is a hard-to-detect 
degenerate object, it is generally difficult to verify such a 
mass-transfer scenario (Pols et al. 1991). The discovery 
of the close degenerate companion to Regulus confirms 
the mass-transfer hypothesis in this case. Considering 
that this a nearby, bright star, it also demonstrates how 
difficult it is to test this scenario in individual cases. In- 
deed, the majority of intermediate-mass stars rotating 
near breakup may have such an unseen companion, sug- 
gesting that further observational scrutiny is warranted. 

2.5. Illustrative Binary Evolution Calculations 

In order to check more quantitatively our proposed 
evolutionary scenario for the formation of the current 
Regulus system, we carried out a number of binary evo- 
lution calculations using a Henyey stellar evolution code. 
All calculations were carried out with an up-to-date, 
standard Henyey-type stellar evolution code (Kippen- 
hahn, Weigert, & Hofmeister 1967), which uses OPAL 
opacities (Rogers & Iglesias 1992) complemented with 
those from Alexander & Ferguson (1994) at low temper- 
atures. We use solar metallicity (Z = 0.02), take a mix- 
ing length of 2 pressure scale heights, and assume 0.25 
pressure scale heights of convective overshooting from 
the core, following the calibration of this parameter by 
Schroder, Pols, & Eggleton (1997) and Pols et al. (1997) 
(for more details see Podsiadlowski, Rappaport, & Pfahl 
2002). 

We ran binary evolution sequences for a substantial 
number of initial mass points lying inside the white re- 
gion of Fig. 2. In most, but not all, cases we adopted a 
value for the angular momentum loss parameter of a = 1. 
The mass retention fraction, (3, is fixed from the initial 
masses and the current binary mass. Since the values 
of (3 typically lie between ~ 0.7 — 1.0, and therefore not 
much mass is lost, our results are not highly sensitive 
to the exact choice for a. We have also verified this by 
running a number of evolutionary models where a = 0.5 
and a = 1.5. The main effect of adopting higher values 
of a is that the initial orbital period can be longer, and 
this tends to push the onset of mass transfer more into 
the regime of early case B mass transfer. 

The model which most closely matches the current-day 
properties of the Regulus system has: M10 = 2.1 Mq, 

4 We note that any of the initial masses specified in Figs. 1 
and 2, transfer more than the requisite amount of mass to spin up 
Regulus. 
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Fig. 3. — Evolutionary parameters as a function of time for the 
most promising model for the current Regulus system. Panel a: 
the evolution of P or b; panel b: the masses of Regulus (long-dashed 
curve), the primary progenitor of the white dwarf (solid curve), 
and the He core of the primary (short-dashed curve); panel c: the 
radius of the primary, the progenitor of the white dwarf; and panel 
d: the mass transfer rate from the primary. 

M 20 = 1.74 Mq, and P orb = 40 hours (1.7 d). This par- 
ticular choice of initial system parameters is shown as a 
filled square symbol within the white region of Fig. [2j 
The first of two sets of results from the binary evolution 
calculation for this model is shown in Fig. [3J The four 
panels show (a) the evolution of P OI b', (b) the progeni- 
tor masses of Regulus (long-dashed curve) and the white 
dwarf progenitor (solid curve); (c) the radius of the pri- 
mary (WD progenitor); and (d) the mass transfer rate, 
M, onto the progenitor of Regulus. Note that as P or b ap- 
proaches 40 days (^1000 hr), the primary has transferred 
most of its envelope mass to the progenitor of Regulus, 
which by then has a mass of ^3.4 Mq. At the same 
time, the primary has developed a ~0.3 Mq degenerate 
He core (short-dashed curve in Fig. 3, panel b) - the pro- 
genitor of the current white dwarf in the system 5 . Fig. [4] 

5 If the progenitor of the white dwarf had been slightly more mas- 
sive and, in particular, the He core mass had exceeded ~ 0.32 Mq 
(which would still be consistent with the present-day measured 
mass function), helium would have been ignited non-degenerately 
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Fig. 4. — Global view of the evolution of the radius of Regulus, 
starting from an epoch just before it starts to accrete matter from 
its companion, the accretion phase (where it reaches its current 
mass of ~3.4 Mq), its subsequent evolution up the giant branch, 
followed by a common envelope (CE) phase. If the system does not 
merge during the CE phase or start mass transfer shortly thereafter 
(see 931, it then becomes a 0.475 Mq core helium burning subdwarf 
and finally a HeCO white dwarf. 

shows the evolution of the radius of Regulus during the 
epoch when it accretes mass from the primary. Note that 
the radius of Regulus, even during the interval when it is 
accreting at ^-10~ 8 Mq yr -1 , never exceeds 4 Rq because 
the mass transfer timescale is substantially longer than 
the Kelvin- Helmholtz timescale. 

2.6. Status of the White Dwarf in the Regulus System 

If the ~0.3 M & companion to Regulus, originally the 
core of the primordial primary in the system, has a mass 
below ^0.31 Mq, then the star will not burn He to 
CO and will cool to become a degenerate dwarf. Af- 
ter a time of several hundred Myr, the maximum time 
allowed since the current Regulus binary has been in ex- 
istence, such a ~0.3 Mq white dwarf would have cooled 
to T eS <15,000 K and a luminosity of <0.02 L G (see, 
e.g., Althaus, Serenelli, & Benvenuto 2001). This is con- 
sistent with the limits on the optical and UV emission 
emission from such a companion in the presence of Regu- 
lus, with a temperature nearly this high and a luminosity 
~10 4 times higher. 

If, on the other hand, the companion to Regulus has 
a mass just slightly in excess of that required for a 
He star to burn He to CO in thermal equilibrium, i.e., 
M ~ 0.32 Mq, then it would still be undergoing nuclear 
burning at the current epoch. The properties of such a 
low-mass He star are: R ~ 0.056 Rq, L ~ 1.4 Lq, and 
T ff ~ 26, 000 K; and the nuclear lifetime is in excess of 1 
Gyr (see, e.g., Hurley, Pols, & Tout 2000). Even though 
such a companion would have a substantially larger lu- 
minosity than its degenerate counterpart, it would still 
not contribute very much light to the Regulus system 
(i.e., < 1%) for wavelengths longward of 1700 A. It even 

in the core (see 82. 6l. 



Evolution of Rcgulus 
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seems unlikely that such a He star would have been de- 
tected in the study of B stars (including Regulus) by 
Morales et al. (2001) using the IUE satellite and EURD 
spectrograph on MINISAT-01. 

2.7. Age of the Regulus System 

One consequence of our scenario is that the age of the 
Regulus system, including any members bound to the 
inner binary, must exceed —900 Myr, the time required 
for the —2.1 Mq progenitor of the white dwarf to evolve. 
Thus, all previous age estimates of — 50 — 200 Myr (see, 
e.g., Gerbaldi et al. 2001) are far too young. The reason 
for these age discrepancies, in essence, is that Regulus 
was rejuvenated by the accretion of nearly 1.7 Mq of 
hydrogen-rich material. 

3. FUTURE EVOLUTION OF THE INNER REGULUS 
BINARY 

3.1. Regulus Ascending the Giant Branch 

In about 100-200 Myr, Regulus will complete its main 
sequence evolution and begin its ascent of the giant 
branch (see Fig. [4] for times after 10 91 yr). The exact 
time to this phase will depend, of course, upon its cur- 
rent core development. The present orbital separation 
of the Regulus binary (at P or b — 40 d) is a ~ 74i? Q , 
and the corresponding Roche-lobe radius of Regulus is 
R L — 44 ± 4.7 R Q . By the time Regulus fills its Roche 
lobe, starts mass transfer to its white dwarf companion, 
and the common-envelope phase commences, it will have 
developed a — 0.48 M Q He core 6 . This is the core that is 
unveiled after the common envelope phase (see the fol- 
lowing section). 

Note, however, that since a 3.4M star will ignite 
helium in its core, its post-CE remnant will continue 
to burn helium. In t he c ase of a successful common- 
envelope ejection (see |3.2| it will appear as an sdB star 
in a very short-period binary for — 10 8 yr (e.g., Maxted 
et al. 2001; Han et al. 2002), a nd ult imately produce a 
hybrid HeCO white dwarf (see ^3.3.2). 



3.2. Common Envelope Phase 

Once mass transfer from the evolved Regulus to its 
white dwarf companion commences, it will quickly be- 
come dynamically unstable (due to the extreme mass 
ratio of the two stars), and a common envelope phase 
will ensue. The relation between the initial orbital sep- 
aration of at ~ 76i?Q, and the post common envelope 
separation, a/ is given by: 



M, 



core , Reg 



M 



Reg 



M 



2M P 



wd 



(4) 



where A is the inverse binding energy of the 
core of Regulus with its envelope, in units of 
-RReg/(GM corC;RGg M cnViRog ), r\ is the fraction of the or- 
bital binding energy that goes into ejecting the common 
envelope, and tli is the size of the Roche lobe of Regulus 

6 We note that whether Regulus fills its Roche lobe during or 
after helium core burning (when it has already developed a CO 
core) depends critically on the amount of convective overshooting. 
For a significantly lower value than our preferred value of 0.25 
pressure scale heights, Regulus would only fill its Roche lobe after 
helium core burning, by which point the hydrogen-exhausted core 
would be considerably more massive (0.7 — 0.8 Mq). 



in units of the initial orbital separation (see, e.g., Web- 
bink 1984; Pfahl, Rappaport, & Podsiadlowski 2003). 
For the parameters of the Regulus-white dwarf binary, 
this expression yields at/aj — 0.0046A?7. The correspond- 
ing post-common envelope orbital period of the binary, 
now consisting of a — 0.5 M Q HeCO dwarf and a 0.3 M Q 
white dwarf, is P or b — 40 (\rf) 3 / 2 minutes. For plausible 
values of Xr/ — 0.1 — 1 (see, e.g., Dewi & Taurus 2000; 
2001; Podsiadlowski, Rappaport, & Han 2003), the post- 
CE orbital periods would likely range between —2 and 40 
minutes, respectively — if both of the compact stars were 
degenerate. Given that the HeCO star is expected to still 
be burning He at the end of the common-envelope phase, 
and that its radius would be no smaller than —0.1 Rq, 
the actual range of post-CE orbital periods is probably 
limited to between —40 and —20 minutes if a merger 
of the cores is to be avoided (see §3.3.1 for more de- 
tails). If the two dwarfs do effectively merge, they would 
produ ce a r apidly rotating giant with unusual properties 
(see, { 3.3.1 1. 



3.3. Post Common-Envelope Phase 
3.3.1. Possible outcomes of the CE spiral-in 

After the common envelope phase, the final orbital 
separation should be a ~ 1/3?7A-Rq, where, as dis- 
cussed above, the product rj\ reflects the efficacy of the 
common-envelope process, and is expected to be in the 
range of 0.1 — 1. What happens to the binary pair consist- 
ing of the 0.3 Mq He white dwarf companion of Regulus 
and the — 0.48 Mq core of Regulus depends critically on 
the uncertain value of 77A. If the 0.3 Mq dwarf is de- 
generate, it will have a radius of 0.017 Rq, and would 
fit well inside its Roche lobe, nearly independent of rj\. 
(If it is still undergoing He burning, its radius would be 
—3 times larger.) The 0.48 M core of Regulus, on the 
other hand, is still burning He to CO in its core, and will 
have a radius of — 0.1 Rq until the He burning phase has 
ended, and the HeCO remnant becomes degenerate (see 
Figs. H] and |6j. This phase lasts for -100 Myr. 

This implies that there is a minimum period of —20 
min for the post-CE system. If the ejection efficiency 
(r/X) is too low, the orbital energy released in the spiral- 
in will not be able to eject the common envelope, and the 
two compact objects, the 0.48 Mq core and the 0.3 Mq 
white dwarf, will merge inside the common envelope. 
This would likely produce a rapidly rotating single gi- 
ant with possibly some very unusual chemical properties 
(perhaps similar to the unusual, rapidly rotating carbon 
star V Hydrae [Kahane et al. 1996]; also see FK Co- 
mae stars [Bopp & Stencel 1981] for related earlier-type 
counterparts) . 

If the common envelope is ejected, the subsequent evo- 
lution depends on the post-CE orbital period, since it 
determines the timescalc on which cither of the two com- 
pact components will achieve contact and start to trans- 
fer mass again. Immediately after a successful common- 
envelope ejection, both components (by definition) will 
underfill their Roche lobes, but the binary orbit will con- 
tinue to decay because of gravitational wave emission. 
The time for orbital decay is 



t = 68 



c.8/3 
orb,CE 



P 



8/3 



Myr 



(5) 



where P Q rb,CE is the orbital period (in hours) immedi- 
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Regulus: Past, Present and Future 



Past: 

stable RLOF 



Present: Q ^Q^p 

He/HeCO WD MS (rapid rotation) 

Future: 

unstable RLOF > dynamical mass transfer 




common-envelope phase 




common-envelope ejection complete merger 

short-period sdB binary with Hc(/CO) WD companion rapidly rotating, single giant (V Hya'.', FK Com?) 



P > 80 min P < 80 min 

AM CVn phase 




degenerate donor He-star donor 

(2 -120 min) (10 -120 min) 



helium nova/SN Ja? 



delayed dynamical merger? 



O sdll binary? 



o 

single sdB star? 

Fig. 5. — Evolutionary scenario: past, present, and future, for the 
Regulus inner binary. Illustrative constituent masses and orbital 
periods corresponding to the different evolutionary phases are: (i) 
Past progenitors - M w ~ 2.3 M ; M 2 o ^ 1.7 Mq; P olb ~ 1 - 15 
days; (ii) Present - M Rog ~ 3.4 Mq; M wd ~ 0.3 Mq; P orb ~ 
40 days; (iii) Future - unstable Roche-lobe overflow; starts with 
system parameters in (ii), but with Regulus having developed ~ 
0.5 Mq He core; (iv) Common Envelope - spiral in of the white 
dwarf into the envelope of Regulus; (v) Post-Common Envelope 
phase — leads to (a) a binary consisting of compact objects in an 
orbit with P or b > 20 minutes (the former core of Regulus now 
appears as an sdB star), or (b) the system merges completely to 
form a rapidly rotating single giant; (vi) Ultracompact Binary - 
depending on the post-CE envelope orbital period, the orbit decays 
due to the emission of gravitational waves to P or b between ~2 and 
~20 min, and the system then appears as an AM CVn system. 



ately after the common envelope, and P or b is the orbital 
period at time t later. The system remains a detached 
binary until either of the two components starts to fill its 
Roche lobe. 

If the immediate post-CE orbital period is < 80 min, 
the 0.48 Mq helium-burning component will start to fill 
its Roche lobe while it is still burning helium in the cen- 
ter, and mass transfer will start from the more massive 
helium-burning star to the 0.3 Mq white dwarf. The sys- 
tem will now become an AM CVn star with a helium-star 
donor (e.g., Nelemans et al. 2 001). The stability of the 



mass transfer is addressed in §3.3.3 

On the other hand 
longer than 



if the post-CE orbital period is 
80 min 7 , the helium-burning star will have 
completed helium-core-burning and become a degener- 



7 Note that such periods are somewhat longer than anticipated 
from our simplistic expression for the ejection of a common enve- 
lope based on energetic arguments (i.e., eq. [4]). 
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Fig. 6. — Evolution of an isolated He star of mass 0.475 Mq. 
The first three panels show the evolution with time of the radius, 
T e g, and luminosity. The mass of the CO core (solid curve; left 
scale) and the central He abundance (dashed curve; right scale) of 
the star are plotted as functions of time in the bottom right panel. 

ate, much smaller HeCO white dwarf (see Figs. [4] and [6| 
before the system becomes semi-detached again. In this 
case, the lighter and physically larger white dwarf will 
fill its Roche lobe first, and mass transfer will take place 
from the lighter white dwarf to the more massive HeCO 
white dwarf. The system will again appear as an AM 
CVn binary (see, e.g., Nelemans et al. 2001), but in this 
case with either a degenerate donor star with Porb — 2 
min and M ~ 0.3 Mq, or with a He-burning star with 
M ~ 0.32 M and P orb ~ 13 min. 

Figure 5 illustrates these different evolutionary paths 
in a "scenario" diagram. 

3.3.2. Evolution of the exposed core of Regulus 

Once the core of Regulus has been exposed after the 
common-envelope phase, it will evolve as an isolated He 
star of mass 0.475 Mq until it or its companion white 
dwarf fill its respective Roche lobe. As mentioned above, 
which star first fills its Roche lobe depends on the or- 
bital period immediately following the common-envelope 
phase. To help understand this, we show in Fig. [6] the 
evolution of an isolated He star of mass 0.475 Mq. The 
first three panels show the evolution with time of the ra- 
dius, T c g, and luminosity. After ^140 Myr, the star has 
completed its He burning, and contracts from ~0.1 Rq 
to its final degenerate radius of ~0.015 Rq. The mass of 
the CO core and the central He abundance of the star are 
plotted as functions of time in the bottom right panel. 
The CO core reaches a final mass of ~0.34 Mq. 

Thus, the critical timescale for the exposed core of Reg- 
ulus to complete its nuclear burning and contract to a 
degenerate state is ~140 Myr. From eq. (5) above we 
see that the original (i.e., post-CE) orbital period would 
have to be in excess of ~80 minutes in order for the 
HeCO star to avoid filling its Roche lobe at P orb ~ 20 
min, and for the system to shrink via gravitational wave 
loses to P or b — 2 min when the He wh ite dwarf would 
first overflow its Roche lobe (see §^3.3.1 and 3.3.3 1. If the 



lower-mass He star is still undergoing He burning (e.g., 
for a mass of ^0.32 Mq) then P orD would be ~13 min. 
rather than 2 min. 

3.3.3. Illustrative calculations of an AM CVn phase 
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Fig. 7. — Illustrative post-CE evolution of the 0.3 Mq white 
dwarf and the HeCO remnant core of Regulus. The post-CE orbital 
period was taken to be 90 minutes. After mass transfer from the 
0.3 Mq white dwarf to the HeCO dwarf commences, the system 
would resemble an AM CVn system with a degenerate donor star. 
The top panel (a) is an HR diagram for the accreting HeCO white 
dwarf. The four panels showing the temporal evolution of the 
system are for (b) the radius of the donor star (He white dwarf, 
solid curve) and the Roche- lobe radius (dashed curve); (c) orbital 
period; (d) constituent masses (white dwarf of initial mass 0.3 Mq, 
solid curve; HeCO star of initial mass 0.475 Mq , dashed curve); and 
(e) mass transfer rate, M. Note that a post-CE orbital period as 
long as 90 minutes is not very probable for the Regulus system. 

In Fig. [7] we show the evolution that results when 
the post-CE orbital period is 90 min. This allows suf- 
ficient time for the HeCO dwarf to complete its He burn- 
ing before gravitational radiation losses bring the sys- 
tem into Roche-lobe contact. The first star to fill its 
Roche lobe and thereby become the donor star is the 
0.3 M & He white dwarf. The initial orbital period at 
the start of mass transfer is ~2min (13 min. if the low- 
mass He star is a bit more massive and is still under- 
going He burning). The mass transfer rate from the 
lower to the higher mass white dwarf is stable. It starts 
at a rate of M ~ 3 x 10~ 6 Mq yr -1 and steadily de- 
clines over the ensuing billion years to < 10~ 12 Mq yr" 1 . 
The details of the evolution of the constituent masses, 
P or b, and M can be seen in Fig. [7] The later parts 
of the binary evolution can, in fact, be computed semi- 
analytically (see, e.g., Rappaport, Verbunt, & Joss 1983; 

Eggleton 2007). We find P orb 
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Mq yr 1 , where ts is the time in 
units of 10 s years from the start of mass transfer. Such a 



Fig. 8. — Illustrative post-CE evolution of the 0.3 Mq white 
dwarf and the HeCO remnant core of Regulus. The post-CE or- 
bital period was taken to be 40 minutes. After mass transfer from 
the 0.48 Mq HeCO dwarf to the 0.3 Mq white dwarf commences, 
the system would resemble an AM CVn system with a helium-star 
donor. Note that P or b at first contact is ~21 min, and then de- 
creases to P or b — 10 min before the orbit starts to expand. The 
various panels, and their description are othewise the same as in 
Fig.0 



system would be observed as an AM CVn star (see, e.g., 
Nelemans 2001). 

In Fig. [8] we show the evolution that results when the 
post-CE orbital period is 40 min. In this case, the or- 
bital decay time is sufficiently short that the post-CE 
~0.48 MQ-helium star does not have sufficient time to 
complete He burning and cool to a degenerate state 
where its radius is smaller than that of the 0.3 Mq He 
white dwarf. In this case, mass transfer starts when the 
helium star is still burning helium in its core (i.e., in 
its hot subdwarf phase), and mass is transferred from 
the hot subdwarf to the 0.3 Mq He white dwarf. Ini- 
tially, mass transfer proceeds on the thermal timescale 
of the helium star due to the fact that the donor is the 
more massive star. However, after an initial phase of 
very high M, the mass ratio quickly inverts, and the 
subsequent mass transfer is driven entirely by the angu- 
lar momentum loss due to gravitational radiation. Note 
that, nuclear burning turns off soon after the helium star 
has lost a significant amount of mass (by this stage, the 
central helium abundance has been reduced from 0.98 to 
0.74, causing the shrinking of the star and the associated 
dip in the mass-transfer rate). The subsequent evolution 
of the binary is similar to that in Fig. [7J except that 
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the now degenerate donor star is helium-depleted and 
oxygen- and carbon-enriched. 

3.3.4. The final fate of the system 

For the binary evolution calculations shown in Figs. [7] 
and [8] any nuclear burning on the accreting star was sup- 
pressed. Unfortunately, the details of the accretion and 
the amount of mass lost from the system in the process of 
accretion and nuclear burning are highly uncertain, and 
at this stage we can only speculate on the final evolution. 
In the AM CVn case with a degenerate donor (Fig. [7]), 
the helium in the accreting HeCO white dwarf will at 
some point reignite in the helium shell under degener- 
ate conditions, producing a helium nova or even a mild 
thermonuclear runaway (possibly producing a ".la super- 
nova" ; see Bildsten et al. 2007) . If helium burning con- 
tinues afterwards, the accretor may swell up significantly 
beyond its Roche lobe, and the system may merge com- 
pletely, producing a single helium-burning hot subdwarf 
(most likely a helium-rich sdO star; Stroer et al. 2007). 

In the case of an AM CVn binary with a non- 
degenerate helium donor (Fig. 181, the accreting white 
dwarf will experience a central helium flash when its mass 
reaches the critical helium-flash mass (~0. 48 Mq). The 
system will probably survive the helium flash, although 
the system may become detached, and the accretor will 
now become a hot subdwarf. 

4. DISCUSSION 

We have discussed quantitatively the interesting past 
and future history of Regulus. The entire evolution- 
ary scenario is summarized in Fig. [5] We have shown 
how the present 40-day orbital period of Regulus and 
its 0.3 Mq white dwarf companion match very well the 
-P or b — M w( j relation predicted by stellar evolution theory. 
The masses of the progenitor stars are inferred to fall in 
the range of M 10 ~ 2.3±O.2M and M 20 ~ 1.7±O.2M . 
Our best numerically computed model is: Mio = 2.1 Mq, 
M 20 = 1.74 Mq, P orb = 40 hours (1.7 d). 



In the future, Regulus will undergo a common enve- 
lope phase wherein the white dwarf will spiral into, and 
eject, the envelope of Regulus. If the envelope of Regulus 
is too tightly bound, then the system may merge, form- 
ing a rapidly rotating single giant with unusual proper- 
ties, possibly related to V Hya and FK Com stars. If, 
on the other hand, the common envelope halts at or- 
bital periods between P or b — 20 and 80 min, the result 
will be a compact binary consisting of a HeCO dwarf of 
mass 0.48 Mq which will start transferring mass to the 
0.3 Mq white dwarf companion (i.e., the current com- 
panion of Regulus). For longer post-CE orbital periods, 
the gravitational wave decay timescale is sufficiently long 
to allow the HeCO star to complete He burning and to 
cool to a largely degenerate state. The system would 
then come into Roche-lobe contact at very short orbital 
periods (P rb — 2 — 13 min). In both of these latter two 
cases, the ensuing mass transfer will likely be stable (at 
least initially) and lead to an AM CVn system whose or- 
bital period, after reaching a minimum, will slowly grow 
to ^20 minutes within ~20 Myr, and to ~40 minutes by 
~200 Myr. Nuclear burning on the accreting dwarf is 
difficult to compute, and we leave this for a future exer- 
cise. Such burning activity could cause the donor to swell 
up, and this, in turn, could cause the compact binary to 
merge. 

Currently we are seeing Regulus as an apparently ordi- 
nary star in the middle of an extraordinary evolutionary 
journey. 
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